%  flit  CC,'"i  AD- A 145  619 


number 


Report  nti;nbeti’i>  MDC  A67  92 


Definition  of  Acceptable  Levels  of 
Mismatch  for  Equivalent  Systems  of 
Augmented  CTOL  Aircraft 


* 


Ls C /:  --'/  'i  '  V  c  n  /?/  C-,-') C 

/>  S'  •  -V.  ’  <<'7  ,  1  7" //"'//  r?  TV 


✓ 

/ 


AfCOOiV/VELL  DOUGLAS 


t  v  /  "7 


$ 


Copy  numbar 


Report  numba&?5'MDC  A6792 


Definition  of  Acceptable  Levels  of 
Mismatch  for  Equivalent  Systems  of 
Augmented  CTOL  Aircraft 

Revision  data  Revision  letter 

Issue  data  19  December  1980  Contract  number 


Prepared 


by  K-W. 


R.  Wood 
Engineer-Technology 
Aerodynamics 


J.  Hodgkin son 

Unit  Chief-Technology 

Aerodynamics 


Approved  by 


R.  r.  Jfehny 
Branch  Chief-Technology 
Aerodynamics 


In  Concurrence--^  ^  % 

C/O’  F.  Mello 


Chief  Technology  Engineer 
Aerodynamics 


EXP 

ELEC 

SEP  1  7 1984 


mcoomwvu.  a/ircivarr  compa/vv 

8»»  5!6,  Saint  Louis,  Missouri  63(66  _  T*l.  (3'4)232-0232 

/ 

IWCOOMJVeu.  DOUGLAS 

'ATVOM 


MDC  A6792 
19  December  1980 


TABLE  OF  CONTENTS 


Section 


Page 


INTRODUCTION 


2.  TWO  FUNDAMENTALS  OF  MANUAL  CONTROL  THEORY 


3.  CRITICAL  ADDED  DYNAMICS 


4.  ENVELOPES  OF  MAXIMUM  UNNOTICEABLE  ADDED  DYNAMICS  ..  . 


5.  MISMATCH  WEIGHTING 


6.  SUMMARY  AND  RECOMMENDATIONS 


7 .  REFERENCES 


APPENDIX  A  DETERMINATION  OF  CTOL  CRITICAL  ADDED 

DYNAMICS . A-l 


APPENDIX  B  CTOL  ENVELOPES  OF  MAXIMUM  UNNOTICEABLE 

ADDED  DYNAMICS . .  .  .  B-l 


APPENDIX  C  MISMATCH  WEIGHTING  .  C-l 


APPENDIX  D  LAHOS  AND  NEAL-SMITH  EQUIVALENT  SYSTEMS  ...  D-l 


LIST  OF  PAGES 


Title  Page 
ii  through  viii 
1  through  22 
A-l  through  A- 2 6 
B-l  through  B-5 
C-l  through  C-5 
D-l  through  D-5 


oouaLAm  conroitArfOM 


MDC  A6792 
19  December  1980 


LIST  OF  ILLUSTRATIONS 


Ficrure 

Page 

1 

Definition  of  Mismatch  . 

.  13 

2 

Attitude  and  Attitude  Rate  to  Stick  Force 
Responses . . . . 

.  14 

3 

Definition  of  the  Crossover  Region  . 

.  14 

4 

Addition  of  Unnoticeable  and  Noticeable 

Dynamics  to  a  Low  Order  System  (Example)  .... 

.  15 

5 

Mismatch  from  Unnoticeable  and  Noticeable 

Added  Dynamics  (Example)  . 

.  16 

6 

VESA  Envelopes  of  Maximum  Unnoticeable  Added 
Dynamics  . 

.  17 

7 

CTOL  Critical  Added  Dynamics  and  Envelopes  .  .  . 

.  18 

8 

CTOL  Envelopes  and  Trams fer  Function  Matches  .  . 

.  19 

9 

CTOL  Envelopes  of  Maximum  Unnoticeable  Added 
Dynamics  (Transfer  Functions)  and  Associated 
Weighting  Factors . 

.  20 

10 

CTOL  and  VESA  Envelopes  of  Maximum  Unnoticeable 
Added  Dynamics  . 

.  22 

A-l 

Block  Diagram  for  Configurations  Simulated  .  .  . 

.  A- 8 

A- 2 

-  A  LAHOS  Added  Dynamics  . 

.  A-9 

A- 3 

-  B  LAHOS  Added  Dynamics  . 

.  A-9 

A- 4 

-  C  LAHOS  Added  Dynamics  . 

.  A-10 

A- 5 

-  2  LAHOS  Added  Dynamics  . 

.  A-10 

A- 6 

-  3  LAHOS  Added  Dynamics  . 

.  A-ll 

A- 7 

-  4  LAHOS  Added  Dynaunics  . . 

.  A-ll 

A- 8 

-  5  LAHOS  Added  Dynamics  ....  . 

.  A-12 

A-9 

-  6  LAHOS  Added  Dynaunics . 

.  A-12 

A-10 

-  7  LAHOS  Added  Dynamics  . 

.  A-13 

A-ll 

-  8  LAHOS  Added  Dynamics  . 

.  A-13 

A-12 

-  9  LAHOS  Added  Dynamics  . 

.  A-14 

iii 

MDC  A6792 
19  December  1980 


Figure 

A-13 

A-14 

A-15 

A-16 

A-17 

A-18 

A-19 

A-20 

A-21 

A-22 

A-23 

A-24 

A-25 

A-26 

A-27 

A-28 

A-29 

A-30 

A-31 

A-32 

A-33 

A-34 


B-l 

B-2 

B-3 

B-4 


-10  LAHOS  Added  Dynamics . A-14 

-11  LAHOS  Added  Dynamics . A-15 

YF-17  LAHOS  Added  Dynamics . A-15 

1A  Neal-Smith  Added  Dynamics . A-16 

6A  Neal-Smith  Added  Dynamics . .  .  A-16 

IB,  2A  Neal-Smith  Added  Dynamics  .  A-17 

6B,  7A  Neal-Smith  Added  Dynamics  .  A-17 

2C  Neal-Smith  Added  Dynamics  .  A-18 

7B  Neal-Smith  Added  Dynamics  .  . . A-18 

7D,  8B  Neal-Smith  Added  Dynamics . A-19 

2E,  3B,  4B,  5B  Neal-Smith  Added  Dynamics  .  ..  .  .  A-19 
6D,  7E,  8C  Neal-Smith  Added  Dynamics . A-20 

IE,  2F,  3C,  4C,  5C  Neal-Smith  Added  Dynamics  .  .  .  A-20 

6E,  7F,  8D  Neal-Smith  Added  Dynamics . A-21 

IF,  2H,  3D,  4D,  5D,  7G  Neal-Smith  Added  Dynamics  .  A-21 

6F ,  7H,  8E  Neal-Smith  Added  Dynamics . A-22 

IG,  2J,  3E,  4E,  5E  Neal-Smith  Added  Dynamics  .  .  .  A-22 

IC,  2B  Neal-Smith  Added  Dynamics  .  A-23 

2G  Neal-Smith  Added  Dynamics  .  A-23 

21  Neal-Smith  Added  Dynamics  .  A-24 

Neal-Smith  and  LAHOS  Critical  Added  Dynamics  .  .  .  A-25 

Definition  of  Maximum  Unnoticeable  and  Minimum 
Noticeable  Added  Dynamics . A-26 

Upper  Gain  Envelope  and  Critical  Added  Dynamics.  .  B-4 

Lower  Gain  Envelope  and  Critical  Added  Dynamics.  .  B-4 

Upper  Phase  Envelope  and  Critical  Added  Dynamics  .  B-5 

Lower  Phase  Envelope  and  Critical  Added  Dynamics  .  B-5 


at coonniu  oouoLM  coKPonxri 


MDC  A6792 
19  December  1980 


Figure 

C-l 

C-2 

C-3 

Table 

A-l 

D-l 

D-2 

D-3 

D-4 


Two  Problems  with  the  Penalty  Method  for  Mismatch 
Weighting . C-3 

Calculation  of  Weighting  Factors  for  an  Arbitrary 
Unnoticeable  Gain  Curve . C-4 

CTOL  Gain  and  Phase  Weighting  Factors . C-5 


Nominal  LAHOS  and  Neal-Smith  Control  System 
Dynamics . A-6 

LAHOS  Equivalent  Systems  -  Standard  Weighting.  .  .  D-2 

LAHOS  Equivalent  Systems  -  Weighting  Factors  .  .  .  D-3 

Neal-Smith  Equivalent  Systems  -  Standard 

Weighting . D-4 

Neal-Smith  Equivalent  Systems  -  Weighting  Factors.  D-5 


las  coesoearfow 


MDC  A6792 
19  December  1980 


CTOL 

PCS 

ps 

L 

a 

MUAD 

n/a 


PR 

S 


v 

'ind 

VESA 


’PCS 


’PS 


5 

SP 

63'  64 


e 


e 


LIST  OP  SYMBOLS  AND  ABBREVIATIONS 
Conventional  Take-Off  and  Landing 
Flight  Control  System 

Longitudinal  stick  force,  positive  aft,  .pounds 

Gain  in  numerator  term  of  longitudinal  stick  deflection 
to  longitudinal  stick  force  transfer  function 

Gain  in  numerator  term  of  pitch  rate  to  longitudinal 
stick  force  bare  airframe  transfer  function 

Dimensional  lift  curve  slope,  sec-1 

Maximum  Unnoticeable  Added  Dynamics 

Change  of  normal  load  factor  with  respect  to  change  in 
angle  of  attack,  g/rad 

Pilot  rating,  Cooper-Harper  scale 

Laplace  operator 

Phugoid  numerator  time  constant,  sec 

Numerator  time  constant  in  pitch  rate  transfer  function; 
approximately  1/L  ,  sec 

Indicated  velocity,  knots 

V/STOL  Equivalent  Systems  Analysis 

Aircraft  elevator  deflection,  positive  trailing  edge 
down ,  rad 

Longitudinal  stick  deflection,  positive  aft,  inches 

Damping  ratio  of  Flight  Control  System  actuator 

Feel  System  damping  ratio 

Phugoid  damping  ratio 

Short  period  damping  ratio 

Damping  ratios  of  control  system  elements 

Pitch  attitude,  deg 

Pitch  rate,  deg/sec 


DOUOL4 

vi 


MDC  A6792 
19  December  1980 


LIST  OF  SYMBOLS  AMD  ABBREVIATIONS  (Continued) 


ti  •  t. 


Lead  or  lag  time  const 2m ts  used  in  control  system 
dynamics ,  sec~l 


Undamped  natural  frequency  of  Flight  Control  System 
actuator ,  rad/sec 


Feel  System  undamped  natural  frequency,  rad/sec 


Phugoid  undamped  natural  frequency,  rad/sec 


Undamped  natural  short  period  frequency,  rad/sec 


U-J  ,  <D, 


Undamped  natural  frequencies  of  control  system 
elements ,  rad/sec 


Subscripts 


Flight  Control  System  actuator 


Feel  System 


Phugoid 


Longitudinal  stick 


Short  Period 


VSLAS  C< 


ifLFi^LifTFgTTnr-r.'  TCT5TT  i**-. 


MDC  A6792 
19  December  1980 


/ 


ABSTRACT 


— ^High  order  pitch  rate  frequency  responses  of  augmented 
aircraft  can  be  broken  down  into  the  response  of  a  low  order 
equivalent  system,  and  the  mismatch  between  the  response  of  the 
actual  high  order  system  and  the  low  order  equivalent.  This 
study  used  variable  stability  NT-33  in-flight  data  to  define 
frequency  response  envelopes  of  acceptable  levels  of  mismatch. 
The  envelopes  were  narrowest  in  a  region  believed  to-  coincide 
with  the  piloted  crossover  region.  The  computer  program  for 
obtaining  the  equivalent  system  was  then  modified  to  weight  the 
match  in  the  crossover  region.  Equivalent  system  summaries  for 
the  LAHOS  and  Neal-Smith  data  sets  are  included. 
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Modem  flight  control  systems  commonly  possess  high  order 
responses.  Flying  qualities  analysis  using  the  requirements  of 
MIL-F-8785B  (Reference  1) ,  which  was  written  in  terms  of  low 
order  classical  aircraft  dynamics,  is  difficult.  One  proposed 
method  of  flying  qualities  analysis  is  to  match  the  high  order 
responses  with  low  order  equivalent  systems,  which  include  an 
equivalent  time  delay.  The  equivalent  system  parameters  are 
then  used  for  assessment  of  the  overall  flying  qualities,  as 
described  in  MIL-F-8785C  (Reference  2). 

A  major  difficulty  with  equivalent  systems  is  that 
acceptable  levels  of  mismatch  between  a  high  order  system  and 
its  low  order  equivalent  have  not  yet  been  defined.  The  mismatch 
has  been  defined  by  MCAIR  as  a  weighted  least  squares  difference 
between  the  gain  (and  phase)  of  the  high  and  low  order  systems 
(Figure  1). 

A  recent  in-flight  simulation  using  the  USAF/Calspan 
variable  stability  NT-33,  sponsored  by  both  the  U.S.  Navy  and 
Air  Force,  examined  the  flying  qualities  of  high  order  systems 
and  their  low  order  equivalents.  A  major  objective  of  this 
experiment,  the  Equivalent  System  Program  (ESP) ,  was  to  define 
the  acceptable  levels  of  mismatch. 

Mismatch  values  of  a  hundred  or  so  proved  unnoticeable  to 
the  pilots  in  the  ESP,  although  previous  work  (References  3-6) 
used  a  mismatch  of  10  as  an  arbitrary  criterion  of  an  acceptable 
fit.  The  criterion  of  10  was  based  on  the  visual  appearance 
of  the  match  when  observed  on  a  Bode  plot;  and  although  it  was 
rooted  in  instinct  rather  than  science,  gave  reasonably  good 
results.  The  insensitivity  of  pilots  to  large  mismatches  in 
the  ESP  requires  an  explanation. 

This  report  offers  a  theory  to  explain  this  insensitivity. 
The  theory  agrees  with  two  fundamental  principles  of  manual 
control  theory.  Following  a  brief  review  of  these  two  fundamen¬ 
tals,  they  are  applied  to  earlier  longitudinal  NT- 3 3  data  on  high 
order  flying  qualities.  Next  the  concepts  of  Unnoticeable  and 
Maximum  Unnoticeable  Added  Dynamics  are  developed.  Finally, 
mcdowwsll  BOUOL4S  conm>*4riow 
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mismatch  weighting  factors  which  ensure  quality  matches  are  cal 
culated  from  envelopes  of  Maximum  Unnoticeable  Added  Dynamics. 
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2.  TWO  FUNDAMENTALS  OF  MANUAL  CONTROL  THEORY 

Studies  such  as  Reference  7  described  the  human  preference 
for  controlled  elements  with  an  attitude  or  attitude  rate  out¬ 
put.  which  is  directly  proportional  to  the  input  command.  The 
longitudinal  short  period  response  of  an  aircraft  can  be  con¬ 
sidered  roughly  an  attitude  rate  system  (Figure  2) .  The  deunping 
and  frequency  limits  of  MIL-F-8785  (References  1  and  2)  can  be 
similarly  considered  to  constrain  deviations  from  the  preferred, 
K/S  response.  Flight  control  systems  can  be  considered  as, 
at  best,  making  the  response  appear  more  K/S-like  (e.g.,  by 
improving  damping  and  frequency) ,  and  at  worst,  introducing 
high  order  effects  which  deviate  from  the  K/S  response  (e.g.,  by 
introducing  equivalent  time  delay) . 

These  studies  have  also  defined  a  crossover  region  of  manual 
control.  This  is  the  frequency  range  in  which  a  K/S-like  response 
is  desired.  For  a  K/S  system,  it  coincides  with  the  unity  ampli¬ 
tude  (0  dB)  crossover  on  an  open  loop  Bode  plot  and  with  the 
closed  loop  bandwidth  as  shown  in  Figure  3.  The  crossover  region 
is  essentially  the  pilot's  frequency  range  of  interest. 
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3.  CRITICAL  ADDED  DYNAMICS 

Typical  high  order  aircraft  dynamics  can  be  represented  by 
some  low  order  or  "base"  system  with  various  lead,  lag,  and  pre¬ 
filter  dynamics  added.  The  addition  of  these  high  order  dynamics 
is  unnoticeable  up  to  a  certain  point,  after  which  the  pilot 
notices  a  change  in  handling  qualities,  evidenced  by  a  Cooper- 
Harper  pilot  rating  differing  from  that  of  the  original  low  order 
system.  Addition  of  high  order  terms  can  improve  a  low  order 
configuration,  because  the  combined  effect  is  that  of  an  improved 
equivalent  system.  Usually,  however,  because  of  the  manual 
control  principles  outlined  in  Section  2,  the  rating  degrades. 

Earlier  NT-33  simulations,  the  Neal  and  Smith  experiment 
(Neal-Smith,  References  8  and  9)  and  the  Landing  Approach  High 
Order  System  experiment  (LAHOS,  Reference  6  and  10) ,  evaluated 
the  effects  of  several  types  of  added  dynamics  on  fighter  handling 
qualities. 

Five  types  of  added  dynamics  (additions)  were  investigated. 

They  were  first  order  lead- lags,  first  order  lags,  second  order 
lag  prefilters,  a  fourth  order  lag  prefilter,  and  a  second  order 
lag  prefilter  -  first  order  lead-lag  combination.  A  range  of 
addition  values  was  tested  within  each  addition  type.  For 
example,  in  the  LAHOS  program,  first  order  lag  addition  time 
constants  (t2)  -1*  .25,  .5,  and  1.0  were  flown.  In  this  present 
study,  the  critical  case  was  defined  for  each  type  of  addition  as 
the  last  set  of  added  dynamics  to  produce  a  high  order  configuration 
that  shows  no  degradation  from  its  base  system's  pilot  rating.  This 
follows  the  approach  of  the  NADC-sponsored  V/STOL  Equivalent  Systems 
Analysis  (VESA,  Reference  11) . 

For  example,  Figure  4  shows  the  frequency  responses  of  a  low 
order  system  (LOS)  along  with  two  high  order  systems  (HOS) 
derived  from  it  by  adding  different  first  order  lead-lags.  The 
configuration  with  the  higher  frequency  lead-lag  addition,  HOS 
#1,  suffers  no  pilot  rating  degradation,  but  the  lower  frequency 
lead-lag  addition  of  HOS  #2  causes  a  definite  degradation.  The 
additions  in  HOS  #1  and  HOS  #2  define  unnoticeable  and  noticeable 
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added  dynamics,  respectively.  In  Figure  5  the  mismatch  from 
Figure  4  is  plotted.  (The  mismatch  is  the  difference  between 
the  LOS  and  the  two  HOS  configurations) .  Based  on  pilot  ratings 
from  Figure  4,  a  HOS  with  mismatch  falling  within  the  HOS  #1 
mismatch  should  have  the  same  pilot  rating  as  the  low  order 
system.  At  some  mismatch  between  HOS  #1  and  HOS  #2,  a  degrada¬ 
tion  in  rating  might  be  expected.  HOS  #1  is  therefore  the 
critical  case  for  the  two  first  order  lead-lag  additions  shown. 

The  process  of  finding  the  critical  cases  was  readily  per¬ 
formed  in  the  VESA  because  the  experiment  was  designed  to 
evaluate  progressive  degrees  of  high  order  contamination  added 
to  low  order  systems  at  various  frequencies.  The  Neal-Smith 
and  LAHOS  programs  were  not  run  with  this  in  mind,  and  therefore 
some  inference  was  needed  when  interpreting  the  data.  Appendix  A 
contains  details  of  the  Neal-Smith  and  LAHOS  critical  cases.  The 
VESA  critical  cases  defined  lateral  V/STOL  envelopes  of  Maximum 
Unnoticeable  Added  Dynamics  (Figure  6) ,  and  similarly  the  Neal- 
Smith  and  LAHOS  critical  cases  define  longitudinal  CTOL  envelopes 
(Figure  7) .  Another  critical  case,  the  LAHOS  phugoid,  was  used  to 
provide  low  frequency  added  dynamics  similar  to  VESA  cases.  In 
Figure  7,  the  frequency  responses  of  the  critical  cases  are  plotted 
on  a  common  Bode  plot  to  obtain  tentative  envelopes  (gain  and 
phase)  of  Maximum  Unnoticeable  Added  Dynamics  (MUAD) . 

As  in  the  VESA,  the  envelopes  defined  the  pilot's  frequency 
range  of  interest  -  the  narrow  region  of  the  mismatch  envelopes. 

In  this  frequency  range  (the  crossover  region  described  previously) , 
the  narrowing  of  the  envelopes  showed  that  pilots  notice  much  lower 
values  of  mismatch  than  at  other  frequencies,  as  would  be  expected. 
Section  4  describes  the  envelopes  in  detail,  and  Section  5  explains 
the  transfer  function  matches  to  the  envelopes. 
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4.  ENVELOPES  OF  MAXIMUM  UNNOTICEABLE  ADDED  DYNAMICS 


The  Neal -Smith  and  LAHOS  critical  added  dynamics  were  found 
and  then  drawn  on  a  common  Bode  plot  as  described  in  the  previous 
section  and  Appendix  A.  However,  the  critical  added  dynamics  were 
all  at  high  frequencies  (i.e.,  at  or  above  the  crossover  frequency) 
unlike  the  VESA  critical  dynamics,  which  included  both  high  and 
low  frequency  additions.  Examination  of  the  VESA  envelopes,  which 
were  wider  at  low  frequencies  than  at  high  ones,  suggests  that 
a  pilot  will  tolerate  relatively  large  gain  or  phase  mis¬ 
matches  at  low  frequencies.  However  a  rational  method  of 
establishing  the  CTOL  envelopes  at  low  frequencies  was  needed. 

The  LAHOS  or  NT-33  phugoid  was  the  sole  contamination  of  the  K/S 
shape  in  this  region,  and  did  not  evoke  adverse  pilot  comments. 

The  phugoid  low  frequency  transfer  function  is: 


phugoid 
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This  phugoid,  the  basic  shape  of  the  VESA  envelopes,  and  con¬ 
siderable  engineering  judgement  were  used  to  shape  the  low 
frequency  portion  of  the  gain  and  phase  envelopes.  The  phugoid 
requirements  of  MIL-F-8785C  (Reference  2)  offered  no  guidance  as 
the  phugoid  frequency  was  not  specified. 

The  envelopes  were  drawn  by  fairing  smooth  curves  either 
through  or  tangent  to  parts  of  the  various  added  dynamics.  The 
four  envelope  curves  (an  upper  and  a  lower  for  gain  and  phase) 
were  tabulated  as  a  function  of  the  frequency  and  then  matched 
by  transfer  functions  using  an  interactive  matching  program 
(see  Appendix  B) .  The  resulting  transfer  functions  define  the 
envelope  boundaries,  are  more  compact  than  tables,  and  can  be 
used  for  any  frequencies  and  frequency  ranges  within  the  limits  of 
0.1  and  100  rad/sec.  These  transfer  functions  are: 
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Upper  Gain  Envelope: 


3.16S  +  31.61S  +  22.79 


S‘  +  27.14  +  1.84 


Lower  Gain  Envelope: 


9.55E-2S*  +  9.92S  +  2.15 


S  +  11. 60S  +  4.95 


Upper  Phase  Envelope: 


68.89S  +  1100. 12S  -  275.22  _.0059S 

- _ -  e 


+  39.94S  +  9.99 


Lower  Phase  Envelope: 


475. 32S  +  184100. S  +  29456.1  .-.0072S 

1 1  it  ■  . -  ■  1  ■  ■  e 


+  11.66S  +  3.89E-2 


These  transfer  fvinctions  provide  good  approximations  to  the 
tentative  envelopes  of  Maximum  Unnotlceable  Added  Dynamics  as 
shown  in  Figure  8,  and  can  be  used  to  examine  the  quality  of 
matches  between  high  order  systems  and  their  equivalents.  A 
mismatch  falling  within  the  envelopes  indicates  that  the  associated 
low  order  equivalent  system  is  satisfactory  as  a  flying  qualities 
tool.  The  envelope  transfer  functions  can  also  be  used  as  a  design 
tool  to  predict  whether  a  set  of  added  dynamics  will  be  noticed  by 
a  pilot. 
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5 .  MISMATCH  WEIGHTING 

The  mismatch  envelopes  allow  determination  of  equivalent 
system  match  quality,  whatever  match  method  has  been  used.  Since 
the  envelopes  imply  that  a  close  match  is  needed  in  the  crossover 
region,  it  makes  sense  to  try  to  ensure  this  close  match  by  some 
means.  Two  possible  methods  of  doing  this  are: 

1)  Penalize  the  match  by  adding  a  large  constant  to  the  mismatch 
function  whenever  the  added  dynamics  fall  outside  the 
envelope;  and 

2)  At  any  frequency  apply  a  weighting  factor  which  is  an  in¬ 
verse  function  of  the  allowable  mismatch  at  that  frequency. 

The  penalty  method  was  rejected  because  of  two  major  faults  which 
would  require  correction  to  ensure  the  desired  match  quality  (see 
Appendix  C) .  It  is  therefore  simpler  to  use  the  second  method 
(weighting  factors)  to  get  good  matches. 

The  mismatch  weighting  factors  can  now  be  calculated  for  any 
frequency  in  a  chosen  range,  because  the  allowable  mismatch  enve¬ 
lopes  are  defined  as  continuous  functions.  In  this  study  the 
chosen  range  was  from  .1  to  100  rad/sec.  The  weighting  factors  are 
based  on  the  allowable  gain  and  phase  mismatch  from  the  Maximum 
Unnoticeable  Added  Dynamics  envelopes,  along  with  the  definition  of 
the  mismatch  function.  Since  the  mismatch  function  at  any  frequency 
was  defined  as  the  square  of  the  difference  between  the  high  and  low 
order  systems'  frequency  responses  at  that  frequency,  the  weighting 
factors  are  similarly  based  on  the  square  of  the  allowable  addi¬ 
tions  at  the  frequency,  as  described  in  Appendix  C. 

The  mismatch  weighting  factors  calculated  from  the  CTOL  gain 
and  phase  MUAD  envelopes  are  shown  in  Figure  9.  The  weighting 
factors  result  from  equating  the  mismatch  along  the  gain  and  phase 
envelope  curves.  The  gain  weighting  factors  at  a  particular 
frequency  were  calculated  by  dividing  the  largest  allowable  gain 
mismatch  by  the  allowable  mismatches  (upper  and  lower)  at  that 
frequency.  The  phase  weighting  factors  were  calculated  similarly. 
Thus  the  weighting  is  largest  when  the  envelopes  are  narrowest. 
Because  the  envelopes  are  not  symmetric  about  the  0  dB  and  0 
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degree  lines,  there  are  separate  upper  and  lower  sets  of  weighting 
factors.  The  choice  of  upper  or  lower  weighting  depends  on  the 
signs  of  the  gain  and  phase  differences  being  weighted.  With 
these  weighting  factors,  good  crossover  region  matches  can  be 
obtained  by  sacrificing  somewhat  the  quality  of  the  match  outside 
the  crossover  region. 


ff 


6.  SUMMARY  AMD  RECOMMENDATIONS 

Using  variable  stability  NT-33  in-flight  data,  it  was 
shown  that  high  order  pitch  rate  frequency  responses  of  augmented 
aircraft  can  be  broken  down  into  two  components: 

1)  The  response  of  a  low  order  equivalent  system  (with  augmented 
frequency,  damping,  and  time  delay  values  which  can  be 
compared  with  MIL-F-8785C  requirements) ;  and 

2)  The  mismatch,  the  difference  between  the  responses'  of  the 
actual  high  order  system  and  the  low  order  equivalent. 

This  study  defined  the  acceptable  levels  of  mismatch  (2) , 
regardless  of  the  method  used  to  obtain  the  equivalent  system  (1) . 
The  method  of  obtaining  the  equivalent  system  was  then  modified 
to  increase  the  likelihood  of  meeting  the  mismatch  criteria. 

The  acceptable  levels  of  mismatch  were  obtained  by  defining 
frequency  response  (Bode)  envelopes  of  Maximum  Unnoticeable  Added 
Dynamics  (MUAD) .  The  likelihood  of  meeting  the  mismatch  criteria 
was  then  increased  by  defining  mismatch  weighting  factors  as  a 
function  of  the  MUAD  envelopes  and  therefore  of  frequency.  The 
envelopes  and  the  weighting  factors  agree  with  the  two  related 
manual  control  theory  concepts  of  the  desirability  of  a  K/S 
response  and  the  importance  of  the  crossover  region;  the  envelopes 
are  significantly  narrower  in  the  crossover  region,  and  similarly 
the  associated  weighting  is  largest  in  the  crossover  region. 

The  envelopes  developed  in  this  study  are  tentative  because 
the  data  were  not  from  programs  specifically  conducted  by 
progressively  adding  high  order  dynamics  to  determine  which  were 
unnoticeable  and  which  were  not.  Some  engineering  judgement  was 
therefore  necessary. 

To  check  or  improve  the  accuracy  of  the  envelopes,  the  effects 
of  various  added  dynamics  could  be  investigated  in  future  simula¬ 
tions,  either  in-flight  or  fixed  base.  The  simulations  could  be 
used  to  expand  or  fill  in  gaps  in  present  configuration  matrices. 
Although  the  work  done  involved  CTOL  longitudinal  dynamics,  a 
similar  approach  to  lateral  and  directional  dynamics  is  warranted. 
Lateral  and  directional  data  for  demanding  tasks  are  lacking. 
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however  the  longitudinal  envelopes  agree  with  general  principles 
of  manual  control  theory  and  are  similar  to  envelopes  developed 
for  lateral  hovering  dynamics  as  Figure  10  shows.  They  are 
therefore  useful  as  tentative  tools  for  lateral  and  directional 
equivalent  system  determination. 

It  should  be  noted  that  the  MUAD  envelopes  allow  for  larger 
mismatches  than  those  normally  seen  in  matching-developed  flight 
control  system  (FCS)  designs.  This  apparently  is  the  .explanation 
for  the  Equivalent  System  Program  (ESP)  results,  in  which  large 
mismatches  were  unnoticeable  to  pilots.  Preliminary  results  also 
indicate  that  the  new  weighting  factors  do  not  alter  equivalent 
system  parameters  significantly  for  those  configurations  which 
were  reasonably  matched  without  the  weighting.  Further  verifi¬ 
cation  of  this  is  required,  however. 
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FIGURE  3  DEFINITION  OF  THE  CROSSOVER  REGION 
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APPENDIX  A 

DETERMINATION  OF  CTOL  CRITICAL  ADDED  DYNAMICS 

Appendix  A  establishes  the  dynamics  used  to  form  the 
envelopes  of  Maximum  Unnoticeable  Added  Dynamics.  The  data 
sets  from  two  conventional  take-off  and  landing  (CTOL) ,  in¬ 
flight  simulations  were  used  to  determine  the  effects  on  handling 
qualities  of  various  types  of  added  dynamics.  Starting  with  the 
lowest  order  or  "base"  configurations  from  the  Neal  and  Smith 
(Neal-Smith)  and  the  Landing  Approach  High  Order  System  (LAHOS) 
programs,  the  critical  cases  for  the  five  types  of  added  dynamics 
or  additions  (first  order  lead-lag,  first  order  lag,  second  order 
lag  prefilter,  fourth  order  lag  prefilter,  and  second  order  lag 
prefilter  -  first  prder  lead-lag  combination)  were  sought.  The 
critical  case  is  defined  as  the  last  set  of  added  dynamics  to 
produce  no  change  in  handling  qualities  or  pilot  rating  with 
respect  to  those  of  the  base  configurations.  The  set  of  all 
critical  cases  is  referred  to  as  the  critical  added  dynamics. 

Table  A-l  defines  all  configurations  in  the  two  data  sets,  and 
Figure  A-l  shows  the  block  diagram  of  the  configurations  simulated. 

Ideally  enough  data  would  be  available  to  allow  easy 
selection  of  a  critical  case  for  each  type  of  addition  (i.e., 
a  finely  varied  matrix  with  progressively  larger  additions  such 
that  several  unnoticeable  additions  would  be  followed  by  an 
addition  that  produces  a  definite  degradation) .  To  account  for 
pilot  rating  scatter,  an  addition  was  considered  unnoticeable  if 
the  worst  degradation  it  caused  was  0.5  (or  less)  on  the  Cooper- 
Harper  scale.  However  the  data  sets  were  of  limited  size  and  there 
were  no  low  frequency  (below  the  crossover  frequency)  added  dyna¬ 
mics.  Sketches  1  and  2  (AGain  Bode  plots)  illustrate  the  procedures 
used  to  find  the  critical  cases  for  the  various  types  of  added 
dynamics. 

The  largest  degradation  for  each  lag  in  Sketch  1  was 
compared  to  the  limit  for  unnoticeability  (0.5)  and  Lags  1  and  2 
were  within  the  limit.  Lag  3  was  not,  so  Lag  2  is  the  critical 
case.  Each  set  of  added  dynamics,  including  the  critical  cases, 
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will  be  referred  to  by  the  configuration  it  degraded  the  most. 
Configuration  A- 2  represents  the  Lag  2  critical  case. 

None  of  the  lags  shown  in  Sketch  2  were  unnoticeable  (within 
the  0.5  limit),  again  based  on  the  worst  degradation  for  each. 
Therefore.no  critical  case  actually  exists  (based  on  the  defini¬ 
tion  given  previously) ,  and  a  replacement  must  be  found.  The 
first  (highest  frequency)  dynamics  of  this  type.  Lag  4  here, 
will  be  used  as  the  critical  case  and  as  before,  will- be  referred 
to  by  the  configuration  with  the  worst  degradation  (Configuration 
C-4  for  the  data  shown  in  Sketch  2) . 

The  degradation  will  serve  as  a  measure  of  how  much  an 
addition  selected  as  a  critical  case  differs  from  its  definition; 
i.e.,  the  smaller  the  degradation,  the  closer  a  "critical  case"  is 
to  being  truly  unnoticeable.  Because  some  "critical  cases" 
actually  are  noticeable,  any  envelopes  developed  from  them  are 
not  conservative;  i.e.,  there  is  a  possibility  that  amismatch 
within  the  envelopes  might  be  associated  with  a  finite  rating 
degradation. 

EFFECTS  OF  LAHOS  ADDED  DYNAMICS 

The  LAHOS  configurations  were  distributed  among  the  five 
types  of  dynamics  added  to  the  base  (-1)  configurations  as 
follows : 

a)  First  order  lead-lag  (-A,  -B,  and  -C) ; 

b)  First  order  lag  (-2,  -3,  -4,  and  -5); 

c)  Second  order  lag  prefilter  (-6,  -7,  -8,  -9,  and  -10); 

d)  Fourth  order  lag  prefilter  (-11) ;  and 

e)  Second  order  lag  prefilter  -  first  order  lead-lag 
combination  (6-)  . 

The  last  type  represented  the  difference  between  the  control 
systems  of  a  modified  YF-1"7  and  the  original  YF-17.  All  con¬ 
figurations  used  the  same  feel  system  and  actuator  dynamics. 

The  LAHOS  critical  cases  for  the  first  order  lead-lag, 
first  order  lag,  and  second  order  lag  prefilter  were  the 
highest  frequency  added  dynamics  of  that  type.  The  critical 
case  was  represented  by  the  configuration  in  this  group  having 
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SKETCH  1 

Frequency 

LAG  1 
LAG  2 
LAG  3 


ADDED 

DYNAMICS 


LAG  1 


LAG  2 


LAG  3 


B-2 :  0. 


B-3:  0.1 


The  data  shown  are  qualitative  and  are  not  related  to  the 
data  sets  actually  used,  except  in  general  format. 


SKETCH  2 
Frequency 


LAG  4 
LAG  5 


ADDED 

DYNAMICS 


[  BASELINE  J  A  [  BASELINE  ]  B  I  [  BASELINE  ]  c 


LAG  4 


LAG  5 


.2 


A-5:  0.8 


B-4 


B-5:  0 


The  data  shown  are  qualitative  and  are  not  related  to  the 
data  sets  actually  used,  except  in  general  format. 
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the  largest  pilot  rating  degradation.  For  the  other  two  types 
only  one  set  of  dynamics  was  used,  and  only  the  configuration 
with  the  largest  degradation  needed  to  be  found. 

The  critical  lead-lag  case  was  represented  by  configuration 
4-C,  which  suffered  a  degradation  of  1.0  on  the  Cooper-Harper 
scale.  The  critical  first  order  ] ag  was  represented  by  2-2  with 
a  degradation  of  2.25.  Although  all  -2  configurations  caused 
degradations,  2-2  caused  the  largest.  The  critical  case  for 
second  order  lag  prefilters  was  represented  by  2-6  (degradation: 

3)  with  the  other  -6  configurations  having  degradations  raging 
from  0.0  to  1.2.  Of  the  fourth  order  lag  prefilter  (-11)  config¬ 
urations,  2-11  and  4-11  had  the  largest  degradation,  6.  The 
second  order  lag  prefilter  -  first  order  lead-lag  combination 
caused  a  degradation  of  8. 

The  frequency  responses  of  the  -C  (first  order  lead-lag), 

-2  (first  order  lag)  and  -6  (second  order  lag  prefilter)  added 
dynamics  were  plotted  on  a  common  Bode  plot  to  begin  the  process 
of  constructing  the  MUAD  envelopes.  The  -11  (fourth  order  lag 
prefilter)  and  yf-17  (second  order  lag  prefilter  -  first  order 
lead-lag  combination)  dynamics  were  not  used  due  to  the  excessive 
degradations  they  caused.  The  LAHOS  added  dynamics  are  shown 
in  Figures  A-2  through  A-15. 

EFFECTS  OF  NEAL-SMITH  ADDED  DYNAMICS 

The  added  dynamics  used  with  the  Neal-Smith  base  configura¬ 
tions  (ID,  2D,  3A,  4A,  5A,  6C,  7C,  and  8A)  consisted  of  seven 
sets  of  first  order  lead-lag  added  dynamics  and  ten  of  first 
order  lag  dynamics.  One  set  of  lead-lag  dynamics  and  two  of 
lag  dynamics  were  also  run  with  a  medium  frequency  Flight  Control 
System  (FCS)  pole  (uiFcg  58  16  rad/sec)  instead  of  the  high  frequency 
FCS  pole  (wpCg  *  63  rad/sec)  used  with  the  other  sets.  The  high 
frequency  pole  had  less  effect  on  the  frequency  response  of  the 
added  dynamics  than  the  medium  frequency  pole,  due  to  the  base 
configurations'  high  frequency  pole  (“FCS  *  75  rad/sec)  nearly 
cancelling  it  out.  The  63  rad/sec  pole  contributed  approximately 
-3  dB  of  gain  and  15  degrees  of  phase  lag  maximum  at  high  frequen¬ 
cies,  while  the  16  rad/sec  pole  added  -26  dB  of  gain  and  100 
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degrees  of  phase  lag. 

Configuration  2C  represented  the  critical  first  order  lead- 
lag  case,  causing  a  degradation  of  0.33.  Lead-lag  configuration 
7fl  was  also  considered  unnoticeable  with  a  degradation  of  0.12. 
Configuration  7D  represented  the  critical  first  order  lag  case 
with  a  degradation  of  2.62. 

The  frequency  responses  of  the  2C  (first  order  lead-lag) 
and  7D  (first  order  lag)  added  dynamics  were  plotted  with  the 
LAHOS  critical  cases  to  aid  in  constructing  the  MUAD  envelopes. 
Although  the  2C  added  dynamics  were  the  critical  lead-lag  case, 
the  unnoticeable  7B  dynamics  allowed  larger  phase  mismatches  at 
high  frequencies.  This  indicated,  at  least  for  the  upper  phase 
curve,  the  envelopes  do  not  "neck  down"  at  high  frequencies.  The 
Neal-Smith  added  dynamics  are  shown  in  Figures  A-16  through  A-32. 

The  critical  cases.  Figure  A-33,  chosen  in  this  appendix  are 
used  in  Appendix  B  to  shape  envelopes  of  Maximum  Unnoticeable 
Added  Dynamics.  The  Maximum  Unnoticeable  Added  Dynamics  is  the 
largest  amount  of  added  dynamics  that  can  be  added  to  a  low  order 
system  without  the  pilot  noticing  a  change  in  handling  qualities. 
This  is  theoretically  quite  different  from  the  Minimum  Noticeable 
Added  Dynamics,  which  is  the  smallest  amount  of  added  dynamics 
that  the  pilot  will  notice.  However  when  examining  real  test  data, 
the  difference  is  hard  to  establish  due  to  the  limited  amount  of 
available  data,  and  the  Minimum  Noticeable  Added  Dynamics  must 
sometimes  be  used  instead  of  Maximum  Unnoticeable  Added  Dynamics. 
Figure  A-34  shows  the  difference  between  the  Maximum  Unnoticeable 
and  Minimum  Noticeable  Added  Dynamics. 
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Data  In  blocka  rafar  cot 


CONFIGURATION /AVERAGE  PR 


CONTROL 

SYSTEM 

DYNAMICS 


SHORT  PERIOD  DYNAMICS 

(nominal) 


ma  * 

n/a  “  4.5  g/rad;  1/t,,  •  .7143  aac-1 


1.0/. 74  I  2.3/. 57  I  2. 2/. 25 


6  2-A  /  5 


1-B  /  5 


2.5  1  3-C 


2-1  /  2 


1-2/5  2-2  /  4.25  3-2  /  7 


1-3  /  9.5  2-3/6 


1-4  /  10  2-4/9 


3. 9/. 54 


5-1  /  6 


1-6  /  5 


1-8  /  8 


2-6/5  3-6  /  6.5 


2-9  /  10 


2-10  /  10 


1-11  /  9  2-11  /  8 


5-3 

/ 

6.2 

5-4 

7 

6 

5-5 

/ 

7 

5-6 

/ 

6 

5-7 

7 

6 

4-10  /  9 


4-11  /_  8  5-11  /  7 


CONFIG. /PR 


6-1  /  10 

YF-17  Original)  |  <.2S+1)(1.1S+1)('  S* 1 2 3  .  2( 

VT7  + 


6-2  '  2  (.5S+1)(.43S+1)(. 065+1) 

YF-17  Modified)  ( . 2S+1) (. 1S+1) (1. 1S+1) 


—  Ill  ■  1  i 

1.9/. 65 

75/. 7 

1.9/. 65 

75/. 7 

(1)  First  number  indicates  base  aircraft  configuration  simulated;  second 
number  or  letter  identifies  control  system  dynamics;  letters  for 
control  system  lead-lag;  numbers  for  lag. 

(2)  Total  configuration  dynamic  model  includes  feel  system  dynamics: 

5 


—  +  2( .6)S  - 

.2  +  26  +  1 


(in/lb) 


(3)  Adapted  from  Reference  10. 

Table  A-l  Nominal  LAHOS  and  Neal-Smlth  Control  System  Dynamics 
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Neal  and  Smith  (Neal— Smith)  Dynamics 


Data  la  blseka  rafar  to: 


comaauxiim/Avatct  n 


SBOBT  retioo  DYMHICS  (Boalnal) 


Vlad  "  250 


s/a  •  18. 3  g/rad;  1/r.  ■  1.2S  aac 


I11BI1JIBU1I 


liiwmii 


1C  /  4.17  2B  /  5. 25 


lUJIflllWi 


NOTES: 


(1)  Number  Indicates  base  aircraft  configuration  simulated;  letter 
identifies  control  system  dynamics. 


(2)  Total  configuration  dynamic  model  Includes  feel  system  dynamics: 

6 

s  .046  , .  „ 

_  .  — - -  (in/lb) 


_  +  K LQIi  +  1 
,2  +  31  +  1 


(3)  Adapted  from  Reference  8. 


Table  A-l  Nominal  LAHOS  and  Neal-Smlth  Control  System  Dynamics  (cont'd) 
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LAHOS  BLOCK  DIAGRAM 
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NEAL-SMITH  BLOCK  DIAGRAM 


Feel  System: 
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Aircraft  Dynamics : 
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Configurations  6-1  and  6-2  had  different  Control  System  dynamics  (up  to  two 
first  order  numerator  roots  and  two  first  order  denominator  roots  in 
addition  to  what  is  already  shown). 


Figure  A-l  Block  Diagram  for  Configurations  Simulated 
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-7  LAHOS  Added  Dynamics 
FIGURE  A-10 
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-8  LAHOS  Added  Dynamics 
FIGURE  A-ll 
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Neal-Smith  Added  Dynamics 
FIGURE  A-30 
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Neal-Smith  Added  Dynamics 
FIGURE  A-31 
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Neal- Smith  Added  Dynamics 
FIGURE  A-32 
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FIGURE  A- 34  Definition  of  Maximum  Unnoticeable  and 

Minimum  Noticeable  Added  Dynamics 
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APPENDIX  B 

CTOL  ENVELOPES  OF  MAXIMUM  UNNOTICEABLE  ADDED  DYNAMICS 

Appendix  B  describes  the  CTOL  envelopes  of  Maximum  Unnotice- 
able  Added  Dynamics  developed  from  the  LAHOS  and  Neal-Smith  criti¬ 
cal  added  dynamics,  the  matching  procedure  used  to  obtain  transfer 
function  matches  for  the  CTOL  envelopes,  and  the  reasoning  behind 
that  procedure. 

Once  all  of  the  critical  cases  have  been  found,  partial 
envelopes  of  Maximum  Unnoticeable  Added  Dynamics  enclosing  them 
may  be  drawn.  The  partial  envelopes  would  enclose  high  frequency 
dynamics  only  and  therefore  be  extremely  narrow  in  the  low  fre¬ 
quency  range.  In  the  VESA  study,  both  high  and  low  frequency 
added  dynamics  were  used  and  this  provided  wide  low  frequency 
envelopes.  The  LAHOS  phugoid  transfer  function: 


S (S  +  1/T.  ) 
91 


S  +  2cp“pS  +  “p 


V 


=  11.77;  u)T 


.196;  Cp  =  .135 


was  added  to  the  Neal-Smith  and  LAHOS  critical  cases  (high  fre¬ 
quency)  to  introduce  low  frequency  effects.  Envelopes  including 
both  high  and  low  frequency  dynamics  were  then  drawn,  Figures 
B-l  through  B-4,  based  on  the  VESA  envelopes.  The  envelopes  were 
input  in  tabular  form  (gain  or  phase  versus  frequency)  to  an 
interactive  frequency  response  matching  program  (NAVFIT)  to  obtain 
transfer  functions  that  match  either  a  gain  or  phase  envelope 
curve  (upper  or  lower)  from  0.1  to  100  rad/sec. 

Using  a  second  over  second  order  (2/2)  system,  the  initial 
match  run  for  each  curve  was  made  with  low  order  time  delay  fixed 
at  zero  and  negative  coefficients  not  allowed.  For  the  same  order 
system,  the  time  delay  and  coefficient  signs  were  allowed  to  vary 
if  a  satisfactory  match  (mismatch  less  than  5)  was  not  obtained 
on  the  initial  run.  If  a  satisfactory  match  still  was  not  obtained, 
the  system  order  was  varied  and  the  procedure  repeated  until  one 
was  obtained . 


WCU.  AfirCM^r 


MDC  A6792 
19  December  1980 


For  the  Upper  Gain  Envelope,  a  2/2  system  with  the  time 
delay  fixed,  positive  coeffients,  and  a  mismatch  of  2.525  was 
chosen.  The  2/2  system  with  time  delay  free  and  negative  coeffi¬ 
cients  allowed  also  had  a  mismatch  of  2.525;  and  a  3/3  system  with 
time  delay  free  and  negative  coefficients  allowed  had  a  mismatch 
of  2.455.  These  systems  were  not  chosen  because  the  initial  2/2 
system  was  simpler  although  the  3/3  system’s  mismatch  was  slightly 
smaller.  For  the  Lower  Gain  Envelope,  2/2  systems  with  time 
delay  fixed  -  positive  coefficients,  and  time  delay  free  -  negative 
coefficients  allowed  were  examined  and  found  to  have  identical 
mismatches  of  1.607.  Using  the  same  reasoning  as  before,  the  time 
delay  fixed  -  positive  coefficient  system  was  chosen. 

The  Lower  Phase  Envelope  was  matched  best  by  a  time  delay 
free  -  positive  coefficient  2/2  system  (mismatch  *  .109) .  A 
similar  system  with  time  delay  fixed  gave  a  mismatch  of  46 .  The 
best  match  obtained  for  the  Upper  Phase  Envelope  had  a' mismatch 
of  96.  This  was  an  unsatisfactory  match  and  the  proposed  envelope 
was  revised.  The  best  match  for  the  revised  envelope  had  a 
better  but  still  unsatisfactory,  mismatch  of  14.6.  After  a  second 
revision,  the  mismatch  for  a  time  delay  free  -  negative  coefficient 
allowed  2/2  system  was  4.8,  a  satisfactory  match. 

The  transfer  functions,  obtained  as  a  result  of  the  matches, 
describe  the  gain  and  phase  envelopes  accurately  and  are  more 
compact  than  tables.  However,  the  major  advantage  of  the  transfer 
functions,  and  the  reason  they  were  found,  is  that  any  frequency 
range  may  be  used  with  them.  If  the  envelope  tables  were  used, 
all  matches  would  have  to  be  performed  at  the  frequencies  used 
in  the  tables.  With  the  transfer  functions,  the  gain  and  phase 
envelopes  can  be  calculated  for  any  frequency  range  and  any 
frequencies  in  that  range.  However,  the  transfer  functions  should 
be  used  between  .1  and  100  rad/sec,  because  the  envelopes  were 
originally  drawn,  tabulated,  and  then  matched  to  transfer  functions 
within  these  limits.  The  transfer  functions  may  be  used  outside 
the  .1  to  100  rad/sec  range  if  the  match  beyond  it  is  checked,  and 
if  necessary  new  transfer  functions  calculated.  The  envelope 
transfer  functions  are: 
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Upper  Gain  Envelope: 

3.16S2  +  31.61S  +  22.79 
S2  +  27.14S  +  1.84 

Lower  Gain  Envelope: 

9.55E-2S2  +  9.92S  +  2.15 
S2  +  11. 60S  +  4.95 


Upper  Phase  Envelope: 

68.89S2  +  1100. 12S  -  275.22  _.0059S 

- 5 - -  e 

+  39.94S  +  9.99 


Lower  Phase  Envelope: 

475. 32S2  +  184100S  +  29456.1  _-.0072S 

- 5 -  e 

S  +  11.66S  +  3.89E-2 
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Figure  B-3  Upper  Phase  Envelope  and  Critical  Added  Dynamics 


Figure  B-4  Lower  Phase  Envelope  and  Critical  Added  Dynamics 
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The  previously  developed  envelopes  of  Maximum  Unnoticeable 
Added  Dynamics  (MUAD)  emphasized  the  desirability  and  importance 
of  a  K  or  K/S  response  in  the  crossover  region.  The  envelopes 
can  not  only  be  used  to  check  match  quality  between  a  high  order 
system  and  its  low  order  equivalent ,  but  also  to  try  to  ensure  a 
high  quality  match.  Two  possible  methods  of  doing  this  are: 

1)  Penalize  the  match  by  adding  a  large  mismatch  constant 
whenever  added  dynamics  or  the  mismatch  between  a 
high  and  a  low  order  system  fall  outside  the  envelope; 
or 

2)  Apply  a  mismatch  weighting  factor  which  is  an  inverse 
function  of  the  allowable  gain  or  phase  mismatch 
(from  the  envelopes)  at  any  frequency. 

The  penalty  method  suffered  from  two  major  faults  that  would 
be  complicated  to  correct.  The  first  fault  was  equal  penalties  for 
any  amount  of  the  frequency  response  outside  the  envelope  (e.g., 
equal  penalties  whether  the  frequency  response  was  1  dB  or  10  dB 
beyond  the  gain  envelope)  as  shown  in.  Figure  Ola.  This  is  easily 
corrected  by  making  the  mismatch  penalty  proportional  to  the 
amount  of  the  frequency  response  outside  the  envelope. 

Even  if  the  first  fault  is  corrected,  equal  amounts  of 
frequency  response  outside  the  envelope  will  be  penalized  equally, 
independent  of  frequency  (e.g.,  a  1  dB  gain  excess  beyond  the 
envelope  at  low  or  high  frequencies  receives  the  same  penalty  as 
a  1  dB  excess  in  the  more  important  crossover  region) ,  as  shown 
in  Figure  C-lb.  This  fault  could  be  corrected  by  a  frequency 
dependent  weighting  factor,  but  this  would  be  similar  to  the  second 
method.  Due  to  the  complications  in  correcting  these  faults,  it 
would  be  simpler  to  use  the  mismatch  weighting  factors  in  the 
first  place. 
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The  mismatch  weighting  factors  were  developed  based  on  the 
definition  of  mismatch  from  Figure  1  and  the  unnoticeable  gain 
and  phase  additions  defined  by  the  envelopes  of  Maximum  Unnotice¬ 
able  Added  Dynamics.  Since  the  envelope  curves  are  the  maximum 
additions  unnoticeable  to  pilots,  and  therefore  their  mismatches 
are  also  the  maximum  unnoticeable,  the  weighting  factors  were 
defined  so  as  to  equate  this  maximum  mismatch.  Calculation  of 
the  weighting  factors  is  illustrated  in  Figure  C-2  for  an  arbi¬ 
trary  gain  response.  Assigning  a  weighting  factor  of  1.0  to  the 
largest  individual  mismatch  resulted  in  weighting  factors  that 
increased  other  mismatches  to  this  value. 

For  each  envelope  there  were  two  sets  of  weighting  factors, 
one  each  for  the  upper  and  lower  curves,  due  to  asymmetric  envelope 
curves.  Figure  C-3  contains  the  weighting  factors  for  the  CTOL 
envelopes  of  Maximum  Unnoticeable  Added  Dynamics. 
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Figure  C-2  Calculation  of  Weighting  Factors  for  an 
Arbitrary  Unnoticeable  Gain  Curve 
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APPENDIX  D 


LAHOS  AND  NEAL- SMITH  EQUIVALENT  SYSTEMS 


This  appendix  contains  four  tables  of  equivalent  systems  for 
the  Neal-Smith  and  LAHOS  configurations.  The  equivalent  systems 
were  obtained  by  matching  the  high  order  pitch  rate  transfer 
functions  of  the  various  configurations  to  low  order  equivalent 
transfer  functions  of  the  form: 


-tS 


Ke<r-+  l>e 

“e 


(D.n2 

SPe 


+  2 


’SP. 


S  +  1 


u> 


SP. 


The  equivalent  systems  were  calculated  using  the  most  recent  ver¬ 
sion  of  the  matching  program,  both  with  and  without  the  weighting 
factors  described  previously.  Tables  D-l  and  D-2  contain  the 
LAHOS  equivalents  (D-l:  Standard  Weighting;  D-2:  Weighting 
Factors)#  and  Tables  D-3  and  D-4  the  Neal-Smith  equivalents  (D-3: 
Standard  Weighting;  D-4:  Weighting  Factors). 
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